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Metals with values of the resistivity and the Hall coefficient much larger than typical ones, e.g., 
of sodium, are called semimetals. We suggest a model for semimetals which takes into account the 
strong Coulomb repulsion of the charge carriers, especially important in transition-metal and rare- 
earth compounds. For that purpose we extend the Hubbard model by coupling one additional orbital 
per site via hybridization to the Hubbard orbitals. We calculate the spectral function, resistivity and 
Hall coefficient of the model using dynamical mean-field theory. Starting from the Mott-insulating 
state, we find a transition to a metal with increasing hybridization strength ("self-doping"). In the 
metallic regime near the transition line to the insulator the model shows semimetallic behavior. 
We compare the calculated temperature dependence of the resistivity and the Hall coefficient with 
the one found experimentally for Yb4As3. The comparison demonstrates that the anomalies in the 
transport properties of Yb4As3 possibly can be assigned to Coulomb interaction effects of the charge 
carriers not captured by standard band structure calculations. 
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I. INTRODUCTION 

Conventional semimetals like the pentavalent elements 
As, Sb and Bi can be well understood using the approxi- 
mation of independent electrons. Having an even number 
of electrons per primitive cell, these solids indeed come 
very close to being insulators. They are not, because 
the occupied band nearest to the Fermi energy over- 
laps slighty with the unoccupied band lowest in energy. 
Therefore a few holes and a few electrons in these two 
bands can contribute to transport leading to a very small 
number of charge carriers n, as seen in the measured 
Hall coefficient i?u = 1/nec of these materials. However, 
the independent electron picture is often not appropiate 
to describe transition-metal and rare-earth compounds. 
Prominent examples are the transition-metal oxides like 
V2O3, MnO or NiO. These oxides are insulators although 
they contain partially filled d bands. Already in 1949, 
Mott identifies the Coulomb interaction of the conduction 
electrons as reason for the insulating behavioiu. Mean- 
while, the term Mott insulator has become a synonym 
for solids with localized charge carriers due to their mu- 
tual Coulomb repulsion. In the early transition-metal 
oxides the movement of the charge carriers is completely 
blocked due to their mutual Colomb repulsion. It seems 
likely that there also exist compounds in which the move- 
ment of the charge carriers is only nearly blocked and 
which exhibit therefore not insulating but semimetallic 
behavior. The aforementioned conventional semimetals 
As, Sb and Bi are solids which come close to being band 
insulators. There might exist transition-metal or rare- 
earth compounds which are semimetallic because they 
come close to being a Mott insulator. In Sec. |^ we will 
discuss a possible example, the semimetal Yb 4 As3. 

In this paper we want to investigate a simple model 
which is able to mimic a semimetal coming close to be- 



ing a Mott insulator. Therefore we are heading for a 
model which can describe the transition from a Mott in- 
sulator to a metal by varying some control parameter. 
The appealing idea is that a metal near a quantum tran- 
sition to an insulator behaves like a semimetal. This 
is so, because the conductivity vanishes for an insula- 
tor. Therefore resistivity and Hall coefficient should be- 
come very large if one approach the insulating state from 
the metallic side. It is important to note, that a non 
vanishing Hall coefficient can be only obtained in a non 
particle-hole-symmetric model. This excludes one pos- 
sible candiate, the well-known Hubbard model with one 
electron per site. Here the control parameter would be 
the ratio between the kinetic and the Coulomb energy of 
the electrons which can be changed by applying pressure. 
A suited model system for our purpose is the self-doped 
Hubbard model introduced in Ref. 0. Here additional 
orbitals are coupled to a half-filled Hubbard model via 
hybridization. The additional electron states serve as a 
generic reservoir of charge. The hybridization drives the 
desired insulator-to-metal transition. 

The self-doped Hubbard model is not based on purely 
academic considerations. It is inspired by the experimen- 
tal observation of charge order in some transition-metal 
and rare-earth compounds like Yb4As3 and NaV205. 
At high temperatures these compounds are mixed- valent 
systems. In Yb4As3 for example the Yb-ions fluctuate 
between a twofold and threefold positive charged state 
leading to a formal valency of Yb 2 25+ . At low tem- 
peratures the charges avoid their mutual Coulomb re- 
pulsion and stay as far apart as possible. In Yb4As3 
this leads to a static charge ordered state in which the 
nearest-neighbor sites of an Yb 3+ -ion are-only occupied 
by Yb 2+ -ions and not by other Yb 3+ -ionsa~Ll. Therefore, 
the Yb 3+ -ion are confined on a subsystem of all possible 
sites, here chains in (lll)-direction. A charge lattice is 
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superposed on the crystal lattice. In YD4AS3 the charge 
order leads to a structural phase transition at T c « 295K. 
The foremost equivalent Yb-sites split up in Yb 2+ -sites 
with a crystal environment of lower symmetry and Yb 3+ - 
site with one of higher symmetry. The superposed charge 
lattice manifest itself in the distorted crystal structure 
below 295K. It is important to note, that charge order- 
ing drive a compound towards a Mott-insulating state 
because the confinement of the charges to a sublattice re- 
duces their kinetic energy. Transition-metal or rare-earth 
compounds whith a charge order transition are therefore 
promising candidates of systems coming close to being 
a Mott insulator. Indeed, Yb4As3 is a semimetaltl. A 
very simplified model have to contain sites preferential 
occupied by the charges representing the sublattice su- 
perposed by charge order. The charges are able to leave 
this sites to gain some kinetic energy. Therefore, the sites 
lower in energy have to be coupled to an empty reservoir 
of charges representing all other accessible but predom- 
inantly unoccupied sites of the crystal. The self-doped 
Hubbard model is a possible realization of this idea. To 
take advantage of sophisticated many-bpdy-, techniques, 
namely the dynamical mean-field theory! 10 !! 11 !, the empty 
reservoir of charges in the self-doped Hubbard model is 
simply represented by unoccupied orbitals hybridizing 
purely local with the ones of the occupied sites. 

The temperature behavior of the resistivity and the 
Hall coefficient is anomalous in Yb4As3. Both quantities 
show a non-monotonous temperature dependence and ex- 
hibit a maximum at a characteristic temperature. The 
characteristic temperatures differ for the two properties. 
Therefore, not only the maximum value of the resistivity 
and the Hall coefficient of the self-doped Hubbard model 
will be of interest but also the temperature dependence of 
the two properties. In this paper we calculate the density 
of states, resistivity and Hall coefficient of the self-doped 
Hubbard model in the vicinity of the metal-to-insulator 
transition using dynamical mean-field theory. In the next 
section we introduce the model and discuss the mapping 
of this multi-band model to an impurity model within 
the dynamical mean-field theory. The impurity model is 
studied numerically by an extension of the non-crossing 
approximation to a two-orbital impurity. We present the 
results in Sec. |TJ, discuss the applicability of the self- 
doped Hubbard model to Yb4As3 
conclude in Sec. [v| 



II. MODEL AND METHOD 

The self-doped Hubbard model is given by 
H = -(A + /i) V/J f, +Ve./J/, -pVcl 



in Sec. IV and finally 
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The model consists of a lattice of / orbitals described 
by the Hubbard Hamiltonian with on-site Coulomb re- 
pulsion U. Locally the / orbitals hybridize with addi- 
tional orbitals called c. The hybridization strength is V. 
For an electron it is energetic favourable to occupy an 
/ orbital. The energy difference between an occupied c 
orbital and an occupied / orbital is given by the charge- 
transfer energy A. The filling is one electron per site, 
i.e., the system is quarter-filled. For V = the self- 
doped model reduces to a half-filled one-band Hubbard 
model. The Coulomb repulsion of the electrons is large, 
U > U c , and the system is a Mott insulator. Here we 
have assumed a finite critical value of the Mott-Hubbard 
transition. For V ^ the mean occupation of / orbitals 
is smaller than one. We may say that for V = the 
electrons order in the sense that they only occupy the / 
subsystem. The charge ordered state is an insulator. A 
finite hybridization hinders the order to be perfect. 

In Ref. H we have shown that for V = 00 the self- 
doped Hubbard model again reduces to a half-filled one. 
However, the effective Coulomb repulsion is reduced to 
[7/4 and the effective hopping to t/2. In a regime of 
U c < U < 2U C we therefore expect an insulator-to- 
metal transition to take place with increasing value of V. 
Unfortunately simple approximation schemes like slave- 
boson mean-field and alloy-analog approximation fail to 
reproduce this quantum transition. In this paper we ap- 
ply the dynamical mean-field theory to the model. 

Though we are aiming at the resistivity and the Hall 
coefficient of the model we first calculate the one-particle 
Green's function G(k,u>). Given the noninteracting 
Green's function, Go(k,oj), of the self-doped Hubbard 
model 
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the self-energies are defined by Dyson's equation 
G{k,u) = [Go L (k,u;)-E(k,oj)' 



(2) 



(3) 



The dynamical mean-field theory assumes a momentum- 
independent self-energy u>) — > T,(oj). In-general the 
self-energy can be considered as a functional of the full 
Green's function G(k, lu). In particular the local approxi- 
mated self-energy is only a functional of the local Green's 
function 
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The functional dependence is generated purely by the in- 
teraction term in the Hamiltonian. It is thus the same for 
an impurity model with the same interactions. Therefore, 
lattice and impurity model have the same self-energy pro- 
vided we identify the Green's function Q(u>) of the im- 
purity model with G(u>). In the case of the self-doped 
Hubbard model the corresponding impurity model reads 
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(i) Gff(u) = £ nm \P? (n,m)| 2 ({X nm ; X mn )) u 

(ii) = ^-/„°° c „dee-' 3e {p„(6)P m ( e + ^)-P i :(e-a;)p m ( e )} 

(iii) ^i mP = E m .r oo de e -^ Pm (e) 

(iv) Pm (e) = -llmP m (e + iO+), P m (0) = de 

(v) P„(z) = 1 (2) 

(vi) S< 2) (z) = -I E mCT delmj(e + iO+) { |Pf (n, m) | 2 /(e)^(* + e) + \P^(m,n)\ 2 f(~e)P m (z - e)} 



TABLE I. Summary of the used equations of the generalized NCA to solve approximatively the impurity problem 
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Note, H c ff just embeds a single unit-cell of the original 
lattice model in an effective medium. Following the dis- 
cussions in Ref. [l^ and 14 it is easy to show that the 
matrix equation 



g{uj) = g(u) 

can be fulfilled by choosing a single function 
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which describes the coupling of the / orbital to the bath. 
The integration in Eq. (||) can be performed analytically 
if we choose a semielliptic density of states 



k 

The self-consistency @j then reduces to 
W 2 



J{u) 



(9) 
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In all our calculation we choose W = 1 as unit of energy. 

What remains is the calculation of the one-particle 
Green's function of the impurity problem (||). Here 
we make use—of an extension of the non-crossing 
approximation!!! (NCA) to the case of more than two 
ionic propagators. Note, that the "impurity" in @ has 
16 eigenstates 
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and not only two as in the case of the U = oo impu- 
rity Anderson model. In Eq. ( pd| ) we have introduced 
the Hubbard operators X nm = |n ){m\. The generalized 
NCA has been applied successfully to the finite-?/ impu- 
rity Anderson modclll the Emery model in the dynami- 
cal mean-field theorytj, and the Andjerson-HubbarcO as 
well as the Kondo-Hubbard modclE-2 in this approxima- 
tion. We just summarize the basic equations in Tab. | 
and refer the reader for further details to the literature, 
e.g. see Refs. ^5| , ^3|Jl4| The coupled integral equations 
(v) and (vi) in Tab. [3 are-solved numerically by intro- 
ducing defect propagatorsEa and making use of the fast 
Fourier transformatiorJl3. The / Green's function of the 
impurity model gives a new function J(lo) via Eq. (Q) 
and the calculations are iterated until self-consistency is 
achieved. 

We are now turning to the calculation of the resistivity 
and the Hall coefficient. In the self-doped Hubbard model 
the electrons can hop only from one / orbital on site i 
to an / orbital on a neighboring site i + 8 a (a = x,y, z 
denotes the direction in space). The current operator of 
the model is therefore given by 

^ =\ at ~ flfi +Sa a) , (^) 

a 

where a is the lattice spacing. In momentum space the 
current operator reads 



(13) 



Given the current operator we can straightforward re- 
peat the calculation of the conductivity in the case of the 
Hubbard model. It is not surprising that we end with the. 
same expressions as in the case of the Hubbard modeO 
just inserting the / spectral density 



fau) (14) 
dw (~E) ftePo&fiM- (15) 



bna 



|e| 3 7r 2 a r°° 



3 




FIG. 1. Spectral density of the / and c orbitals in the 
self-doped Hubbard model for U — 3, A = 2 and different 
values of the hybridization V. For V > 2.7 the system is 
metallic and a sharp peak at the chemical potential is present 

Here f(u>) = [exp(/3w) + 1] 1 is the Fermi function. The 
derivation involves two steps. First, we consider the non- 
interacting case like in Ref. [HI. Second, we mark that 
vertex corrections in the linear response diagrams van- 
ish. Here the fact enters that the one-particle self-energy 
is momentum independent. In addition the special k- 
dependence of the free propagators via eg and of the ver- 
tices via (see Eq. (|l3|)) is needed.0 Because the ver- 
tex corrections vanish the free propagators in the expres- 
sion for the noninteracting case can simply be replaced 
by full ones. 

Given the conductivity a xx and the Hall conductivity 
a? R Z we know the resistivity p xx — l/a xx and the Hall 
coefficient i?H = cr^" 2 ' / (a xx ) 2 . The units of the two trans- 
port properties are given by [p xx ] ~ [ha/c 2 ] « 0.1mi7cm 
and [Rh] « [a 3 /ec] « ilO _5 cm 3 /C, respectively. Note, 
the conductivities are given by the one-particle spectral 
function only 

p f {e,uj) = --lmGff{e,w + iO + ) . (16) 

In particular, if pf(e,p) vanishes for all energies e, i.e., 
Pf(p) = 0, than also the conductivity vanishes for low 
temperatures. In this case the system is an insulator. 

III. RESULTS 

When discussing the results we will frequently refer 
to experimental findings for the semimetal Yb4As3. In 
the next Section we will discuss the applicability of the 
self-doped Hubbard model to this rare-earth compound. 




FIG. 2. Spectral density of the / and c orbitals in the 
self-doped Hubbard model for U = 5, A = 2 and different 
values of the hybridization V. The system stays insulating 
independent of the value of V. 

Heading for the metal-insulator diagram of the self- 
doped Hubbard model we first concentrate on the value of 
the / spectral function at the chemical potential p / (/i) for 
different values of the hybridization. In Fig. |l] we display 
the evolution of the spectral function of the model with 
increasing V and fixed values for U and A. The / and c 
spectral functions are obtained by the treatment outlined 
in Sec. Q. Within our approach the critical Coulomb 
repulsion of the Mott-Hubbard transition in the one-band 
Hubbard model is t/ c — 1.77(5). The chosen value of U = 
3 fullfills U c < U < 2U C . As expected for this value of 
U the model undergoes an insulator-to-metal transition. 
For V < 3 we obtain pf(p) = and the system is an 
insulator. For V > 3 we find a finite value Pf{p) ^ and 
the model is a metal. Increasing the resolution in V we 
obtain the value V c ~ 2.7 for the critical hybridization 
of the transition. Next, we consider the case of a very 
large Coulomb repulsion {7 = 5, i.e., U > 2U C . In Fig. 
U we demonstrate that the gap between the filled band 
and the lowest empty band never closes. The density of 
states at the chemical potential vanishes for all values of 
the hybridization. This shows clearly that the insulator- 
to-metal transition is restricted to a parameter region 
U c < U < 2U C . In Fig. |^ we display the resulting metal- 
insulator diagram of the self-doped Hubbard model in 
the (U, y)-plane for two different values of the charge 
transfer energy. 

Qualitatively a phase diagram of this form was al- 
ready predicted in Ref. |[ Indeed, the bands in the spec- 
tral function of the lattice model evolve like the tran- 
sition energies of the single two-orbital impurity (H cc u 
in Eq. (H)). Therefore no sophisticated calculations are 
required to determine roughly the position of the cen- 




FIG. 3. Metal-insulator diagram of the self-doped Hubbard 
model for A = 2 (full circles) and A = 1.5 (open circles). 
The lines are guide to the eye and separate the metallic from 
the insulating regime. Transitions between the two states 
are limited to a region between U c = 1.775 and 2U C = 3.55 
marked by horizontal dashed lines. 

ters of bands. We denote the possible transitions of the 
single two-orbital impurity by (f^c ) — > (/™c m ). The 
correspondence of the respective transition energies and 
the band energies allows us to classify the occupied band 
below the chemical potential as (/°c°)-band and the un- 
occupied bands above \x as (/ 1 c 1 ) s -, (/ 1 c 1 ) t - and (/ 2 c )- 
bands (here s and t refers to singlet and triplet). Lat- 
tice effects show up in the ratio between / and c weight 
constituting a single band. For example, the occupation 
number of c orbitals is larger in the lattice model than 
in the two-orbital impurity marking a larger c weight in 
the occupied (/°c°)-band. This is a simple consequence 
of quantum fluctuations between different impurity con- 
figurations in the ground state of the lattice model. 

Nevertheless, the most prominent feature of the lattice 
model is the appearance of a sharp Kondo-likc resonance 
close to the chemical potential. The temperature depen- 
dence of the resonance is shown in Fig. ||. The peak arises 
below a characteristic temperature To. Regarding the hy- 
bridization dependence of the resonance we state that it 
becomes narrower and shows up at a lower temperature 
Tq(V) when we approach the insulating state from the 
metallic side, i.e., when we regard the limit V — > . 
Note, that for temperatures just above To /0/(/i) is a 
monotonously decreasing function in this limit. There- 
fore the observed behavior of To is very similar to the 
well-known dependence of the Kondo temperature on the 
density of states Tk ~ exp (—b/p(fi)). The usual Kondo 
effect cap be interpreted in terms of the resonant-level 
model.E3 The conduction electrons with energies close to 
/i are virtually bound to localized electrons due to res- 
onance scattering. In a similar way one may assign the 



FIG. 4. Temperature dependence of the / spectral den- 
sity of the self-doped Hubbard model with parameters 
(U,A,V) = (1.8,2,0.6). The inset shows the enlarged fre- 
quency region close to the chemical potential. Below the tem- 
perature To a sharp peak arises at the chemical potential. 

resonance seen in our calculations to a band-Kondo ef- 
fect where the conduction electrons are virtually bound 
not to localized but to moving electrons. This "band- 
Kondo effect" is a typical feature of strongly correlated 
metals treated within the dynamical mean-field theory. 
It occurs also in the one-band Hubbard model. It is still 
a matter of dispute if the effect is physical or a short- 
coming of the mapping of the lattice model to an im- 
purity model as done within the dynamical mean-field 
theory. But at least in the case of the Hubbard model 
there are indications that the effect can be seen also by 
other methods .EiJE3 

The strong temperature dependence of the spectral 
density close to the chemical potential causes a strong 
temperature dependence of the resistivity. In Fig. |5| we 
show this dependence as obtained by the treatment out- 
lined in Sec. ||. First we focus on the magnitude of the 
resistivity. The resistivity increases when we decrease 
the hybridization, i.e., when we approach the insulat- 
ing regime. In fact we expect a diverging resistivity in 
the limit V —> because the conductivity vanishes in 
an insulator. This is the key idea. In the vicinity of a 
quantum transition to an insulator we can obtain large 
values for the resistivity and the Hall coefficient implying 
semimetallic behavior. In Fig. || the order of magnitude 
of the resistivity is lmf2cm, i.e., the same as observed in 
Yb4As3. Note, that large values of p and Rn do not imply 
that the occupation number of the c orbitals rf = 1 — n* 
have to be small. Large values are possible for every n c 
in striking contrast to a simple Drude picture when as- 
suming the charge carriers in the system are given by 
the missing electrons in the / subsystem. We stress this 
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FIG. 5. Temperature dependence of the specific resistivity 
in the self-doped Hubbard model. The unit of the the resis- 
tivity is O.lmficm. Magnitude and position of the maximum 
shift with the value of the hybridization V. 

point because in YD4AS3 there is experimental evidenceS 
that the number of missing holes in the (lll)-chains is 
not identical with the low carrier concentration obtained 
from the large Hall coefficient. 

As function of temperature the resistivity shows a max- 
imum at a characteristic temperature Tq. This temper- 
ature is identical with the one discussed above, i.e., the 
temperature where the sharp resonance arises in the spec- 
tral function. As seen in Fig. || and as discussed before 
To is larger for parameter values deeper in the metallic 
regime of the self-doped Hubbard system. Yb4As3 can be 
tuned to a more insulating or metallic state (indkjated by 
the value of the resistivity, bv applying pressures or sub- 
stituting P or Sb for As.EjO Indeed, the experimentally 
observed characteristic temperature shows the expected 
behavior. In samples with a lower (more metallic) resis- 
tivity To becomes larger as compared to samples with a 
higher resistivity. Note, that the resonance appears also 
in the c spectral function (see Fig. [l] or Fig. ^) . Con- 
cerning Yb4As3 we may interpret the c orbitals as As 
p band with zero bandwidth. However, in Yb4As3 the p 
band is broad and especially the holes in this band should 
contribute to transported. We may conjecture that also 
the p holes reveal the described low temperature scale 
T (see also the discussion in the next Section). This is 
known for the case of a periodic Anderson model where 
also a non-monotonous, temperature dependence of the 
resistivity is obtained.E3 

The enhanced / spectral function close to the chemical 
potential indicates heavy masses of the charge carriers at 
low temperatures. Unfortunately, our approach does not. 
allow us to perform the limit T — > 0. It is wcll-knownc3 
that the non-crossing approximation underestimates the 



FIG. 6. Temperature dependence of the Hall coefficient in 
the self-doped Hubbard model. The unit of the Hall coefficient 
is ilO _5 cm 3 /C. Magnitude and position of the minimum 
shift with the value of the hybridization V. 

absolute value of the imaginary part of the self-energy of 
the single-impurity Anderson model. This failure will be 
enhanced in the self-consistent adjustment of the hosting 
bath of conduction electrons leading to an unphysical 
change of sign of the self-energy at low temperatures. E£l 
Therefore we cannot proof the relation p(T) = AT 2 and 
we cannot calculate the coefficient A. However, note that 
in Fig. |5| the onset of a T 2 -behavior can be seen at least 
for the case of V = 0.6. 

We now turn to the temperature dependence of the 
Hall coefficient as displayed in Fig. ^. For discussion we 
want to compare qualitatively the calculated Hall coeffi- 
cient with the one measured in Yb4As 3 . First of all the 
Hall coefficient of the self-doped Hubbard model is neg- 
ative whereas it is positive in Yb4As3. From our point 
of view this has to be expected. In the model electrons 
order in the sense that they prefer to occupy the / sub- 
system. In Yb4As3 however holes in the 4/-shells of the 
Yb-ions are predominantly confined to the (lll)-chains. 
So there is a different sign for the charge carriers which 
should show up in the Hall coefficient (see also the mod- 
ified model in the next Section). Second, the magnitude 
of the Hall coefficient seems to be wrong. In Fig. || it 
is still five orders of magnitude too small. This rather 
small value of the Hall coefficient is caused by the still 
very symmetric spectral function close to the chemical 
potential (see Fig. ||). In this respect the chosen values 
of the hybridization (about 25% of the unperturbed / 
bandwidth 2W) are too large. Nevertheless, as we have 
stressed before the value of the Hall coefficient can be 
tuned by choosing parameter values of the model closer 
to the metal-to-insulator transition. Close to the transi- 
tion arbitrary large values for the Hall coefficient may be 
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obtained. 

Keeping this in mind we find indeed an anomalous tem- 
perature dependence of the Hall coefficient as in exper- 
iment with a characteristic temperature TMax- In the 
case of V = 0.6 in Fig. || this characteristic tempera- 
ture is smaller than the one found in the resistivity (see 
Fig. ||). The relation Tq/Tm^k ~ 2.13 is close to the ex- 
perimental one Xo/Xiviax « 1.75. As we will point out 
in the next Section, the self-doped Hubbard model is 
not a model simplifying the band structure of Yb4As3 
close to the chemical potential. Therefore, we only want 
to discuss possible implications of the Coulomb interac- 
tion of the charge carriers on the transport properties. 
If we would e.g. identify the temperature To = 0.0425 
with the experimental one of 140K we obtain for the 
/ bandwidth 2W « 0.6eV. The LSDA+t/ approach 
which will be reviewed in the next Section gives a value 
of 2W = 0.007eV. But even qualitatively there is still 
an essential difference. Surprisingly XJyiax decreases with 
increasing value of the hybridization. This behavior of 
^Max is just opposite to the one of T . In experiment 
both T and Xfviax behave the same. Especially the ratio 
p(T)/Rn(T), i.e., the inverse Hall mobility of the charge 
carriers, is independent of the applied pressure. It is not 
clear if our finding for the behavior of TMax in the model 
depends on the parameter regime and possibly change 
closer to the metal-to-insulator transition or not. 



IV. THE SEMIMETAL YB 4 AS 3 

In the preceding Section we have presented the 
semimetal Yb4As3 as possible candidate for a compound 
which can be described by the self-doped Hubbard model. 
In this Section we will describe the properties of this com- 
pound in more detail. We want to argue that it is indeed 
justified to compare the transport properties of the self- 
doped Hubbard model and of Yb4As3 as we have already 
done. 

Yb4As3 is an example of a low-carrier heavy-fermion 
systemEI. Below 100K a linear specific heat is found with a 
large coefficient 7 ss 200 mJ/(molK 2 ). The Sommcrfcld- 
Wilson ratio is of order unity. The large Hall coefficient at 
low temperatures indicates an extremely small concentra- 
tion of 6 = 0.3 x 10~ 3 positive charge carriers/Yb-atom. 
The large Hall coefficient and the large residual resistiv- 
ity of 1 milcm at low temperatures classify Yb4As3 as 
a semimetal. In Yb4As3, all Yb-atoms are aligned on 
four families of chains. Assuming trivalent As the ratio 
of Yb-ions is Yb 3+ : Yb 2 + = 1 : 3. A Yb 3+ -ion has one 
hole in the / shell (4/ 13 ) and shows a magnetic moment 
due to Hund's rule coupling. In a sequence of papersETTJ 
Schmidt, Thalmeier and Fulde suggests that at low tem- 
peratures the Yb 3+ -ions are confined to one of the four 
chain systems, e.g., parallel (lll)-chains. The low energy 
excitations of these quasi-one-dimensional spin chains ex- 
plain the specific heat including the large 7-coefficient. 



However, due to the strong Coulomb repulsion of the 
holes and the large distance of the Yb 3+ -sites within one 
chain a perfect ordered state would imply an insulator. 
Because Yb4As3 is a semimetal a small fraction of holes 
have to be redistributed either on the As-atoms or on 
the other three chain systems, i.e., the (lll)-chains have 
to be "self-doped" . It is tempting to identify this small 
fraction of holes with the low concentration of charge car- 
riers seen in the Hall coefficient. The charge orderingppi 
the Yb 3+ -ions is meanwhile experimentally confirmed-Lru 
However, the number of missing holes in the (lll)-chains 
cannot be exactly measured. Note, that the polarized 
neutron diffraction experiment suggest that this number 
is not small but of the order of several percent. 

During the completion of our work, the electronic 
structure of Yb4As3 has been investigated using en- 
ergy band- calculations within the so-called LSDA+£/ 
approached. The calculations take into account the pos- 
sibility of hole occupation on the Yb-sites of the (111)- 
chains and on the As-sites but not on the remaining three 
chains system. In result a very narrow Yb 3+ 4/ hole band 
is obtained close to the top of a broad As p band. The 
Fermi energy is pinned to the bottom of the 4/ hole band. 
From this band structure indeed a very small occupation 
number of holes on the As-sites results comparable to the 
small carrier number estimated from the Hall coefficient. 
This finding is independent from the chosen value of U 
over a wide range of several eV. One should mention that 
the possible hole distributions in the used approach are 
too restrictive to deal with valence fluctuating systems 
like Yb4Bi 3 and Yb4Sb 3 which do not show a charge or- 
dering transition. Only the p holes should contribute 
to transport properties because there is no direct //- 
hopping. This explain correctly the sign and magnitude 
of the Hall coefficient. 

Some questions remain intriguing. In agreement with 
Fermi-liquid behavior the resistivity is found to be of 
the form p(T) = po + AT 2 with a large coefficient 
A w 0.84/iilcm/K 2 . It is hard to understand which scat- 
tering mechanism of the p holes can lead to such high 
effective masses in transport. Not understood at all is 
the temperature dependence of the resistivity and the 
Hall coefficient of Yb4As3. Both quantities show a non- 
monotonous temperature dependence and exhibit a max- 
imum at a characteristic temperature. The characteristic 
temperatures differ for the two properties. 

Two comments seem to be appropiate. First of all, 
the band-structure results question strongly the validity 
of the self-doped Hubbard model for the compound un- 
der consideration Yb4As3. In this compound the trans- 
port seems to be purely due to the few charge carriers in 
the nearly empty charge reservoir, e.g., the holes on the 
As-sites. In the self-doped Hubbard model however the 
electrons in the c-states do not contribute to the trans- 
port properties of the model at all. Moreover, the per- 
fect charge ordered state, e.g., all holes confined to the 
(lll)-chains, is a Mott insulator of a very extreme type 
because the kinetic energy of the holes on the Yb-sites 
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FIG. 7. Spectral density of the / and c orbitals in the 
self-doped Hubbard model for U — 3, A = 2 and V = 4. In 
comparison with Fig. ji] the enlarged frequency region close 
to the chemical potential is shown. 

is nearly zero. Therefore, a suitable model Hamiltonian 
is still of the form of Eq. (|l|) but with two important 
changes. The filling should be three electrons per site 
and the Coulomb repulsion should be suffered by the o 
electrons. In this modified version the /-subsystem (and 
not c) represents the As sites with nearly no holes and the 
c-subsystem represents the Yb 3+ -ions occupied by nearly 
one hole per site. In a future work we will investigate this 
model. Nevertheless, we believe that the modified model 
behaves quite similar to the self-doped Hubbard model. 
Following the arguments outlined in Ref. an insulator- 
to-metal transition as function of the hybridization have 
to be expected for both, the self-doped Hubbard model 
and the modified model. In Fig. [?| we show the / and 
c spectral functions of the self-doped Hubbard model in 
the metallic regime. A sharp peak close to the chemi- 
cal potential is seen in both, the / and the c spectral 
function. The interaction induced many-body effects of 
the correlated subsystem are carried over to the system 
of the uncorrelated orbitals. Therefore, the similarity of 
the self-doped Hubbard model with the proper model of 
Yb4As3 holds despite the fact that the conductivity of 
the modified model is given by the spectral function of 
the uncorrelated orbitals. In conclusion, it is justified to 
compare here some findings for the self-doped Hubbard 
model with transport measurements of Yb4As3. 

We can argue in a slightly different way. In the charge 
ordered phase of Yb4As3 the Coulomb repulsion con- 
strains the holes in the 4/-shell of the Yb-ions to be 
localised. The long range part of the repulsion restricts 
the holes to chains in (lll)-direction. The on-site part 
hinders the holes to move along a single chain. This im- 
plies that the 4/-spectral density is split into bands sepa- 
rated by gaps which we may call charge order and Mott- 



Hubbard gaps, respectively. The charge order gap will be 
of the order of the transition temperature T c tzs 295K and 
is much smaller than the Mott-Hubbard gap of several 
eV. Therefore, the relevant gap for electron excitations 
is the charge order gap and not the Mott-Hubbard gap. 
However, the charge order gap in the 4/-spectral density 
is a true many-particle effect and cannot be explained by 
the change of the unit cell between the charge disorderd 
and the charge ordered phase. Note, the size of the unit 
cell is unchanged in the charge order transition. The cell 
is only distorted, changing from cubic to trigonal. The 
distortion of the unit cell has only little influence on the 
electronic structure of Yb4As3, as shown in Ref. ^6|. We 
believe that the charge order gap resembles much more 
a Mott-Hubbard gap than it resembles a gap in an usual 
band insulator. The physics of the charge ordering is not 
included in the self-doped Hubbard model. But what is 
considered in the used model is the effect of hybridiza- 
tion between a correlated (Hubbard like) band and an 
uncorrelated band. 

The second comment concerns the temperature de- 
pendence of the resistivity and the Hall coefficient of 
Yb4As3. On the first sight it resembles the one found in 
other heavy- fermion compounds like UPt3. In UPt3 the 
anomalous Hall conductivity arises from skew-scattering 
of the charge carrier aEI. As a consequence, the anoma- 
lous Hall coefficient should be proportional to the square 
of the resistivity for temperatures smaller than Tjviax, 
the temperature for which a maximum in the Hall co- 
efficient is observed. However, in Yb4As3 the relation 
i?H ~ p 2 is not fullfilled in the temperature regime from 
80K down to 4K. Moreover, the zero-temperature value 
of the Hall coefficient in Yb4As3 is at least three orders 
of magnitude larger than in other heavy-fermion com- 
pounds with a positive Hall coefficient. We therefore ex- 
clude skew-scattering as possible mechanism to explain 
the experimetal data. 

V. CONCLUSION 

In conclusion, we extended the Hubbard model by cou- 
pling one additional orbital per site via hybridization to 
the Hubbard orbitals. The Coulomb repulsion is sup- 
posed to be large, U > U c , where U c denotes the critical 
Coulomb repulsion for the Mott-Hubbard transition in 
the one-band Hubbard model. We calculated the spec- 
tral function, resistivity and Hall coefficient of the "self- 
doped" Hubbard model using dynamical mean-field the- 
ory To this end the lattice model is mapped onto an 
impurity model in which a unit cell of the lattice is em- 
bedded self-consistcntly in a bath of free electrons. The 
impurity model is studied numerically by an extension 
of the non-crossing approximation to a two-orbital im- 
purity. 

The self-doped Hubbard model is an insulator only in 
a restricted parameter regime. The hybridization with 
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the added orbitals drives an insulator-to-metal transi- 
tion, provided the Coulomb repulsion fullfills the con- 
straint U c < U < 2U C . This is a correlation driven 
metal-insulator transition in a non particle-hole symmet- 
ric case. In the vicinity of the transition to an insulator 
the resistivity and the Hall coefficient become very large. 
Therefore the model serves as a model for semimetallic 
behavior in systems where the Coulomb interactions of 
the charge carriers nearly block their movement. The 
number of missing electrons in the Hubbard subsystem, 
i.e., the occupation number of the additional orbitals, can 
not be interpreted as number of charge carriers seen in 
the transport properties. A simple Drudc analysis fails 
for this type of mechanism leading to semimetallic be- 
havior. 

In the semimetallic regime the resistivity and the Hall 
coefficient show an unusual temperature dependence. 
Both, p(T) and |.Rh(T)|, exhibit a maximum at a char- 
acteristic temperature. The characteristic temperatures 
for the resistivity To differs from the one of the Hall coef- 
ficient Tjviax- Depending on the parameter regime either 
To < Tviax or T > Tviax is possible. The characteristic 
temperatures change considerably when varying the hy- 
bridization and tuning the system into a more insulating 
or metallic state. At To a resonance arises in the spectral 
function close to the chemical potential. This resonance 
can be assigned to a Kondo-likc effect for band electrons. 
The enhanced spectral weight at the chemical potential 
indicates heavy masses for the charge carriers. 

We compared our findings with measurements of the 
resistivity and the Hall coefficient for the semimetal 
Yb4As3 including pressure and substitution experiments. 
Although the self-doped Hubbard model is not a sim- 
plified model extracted from band structure calculations 
we found surprising similarities. We take this as indica- 
tion that the anomalous transport properties of Yb4As3 
can be indeed assigned to the Coulomb interaction of the 
charge carriers which is not treated adequate in stan- 
dard band structure calculations. This should be true 
even when the transport properties are dominated by As 
p holes. At least in the model the interaction induced 
many-body effects of the Hubbard subsystem are carried 
over to the system of the additional orbitals. Neverthe- 
less, to obtain final conclusive results the band structure 
of Yb4As3 have to be used as input for a simple model 
which still can be treated within the dynamical mean- 
field theory. This is left for future work. 
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